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We calculate contained and upward muon flux and contained shower event rates from neutrino 
interactions, when neutrinos are produced from annihilation of the dark matter in the Galactic 
Center. We consider model-independent direct neutrino production and secondary neutrino pro- 
duction from the decay of taus, W bosons and bottom quarks produced in the annihilation of dark 
matter. We illustrate how muon flux from dark matter annihilation has a very different shape than 
the muon flux from atmospheric neutrinos. We also discuss the dependence of the muon fluxes on 
the dark matter density profile and on the dark matter mass and of the total muon rates on the 
detector threshold. We consider both the upward muon flux, when muons are created in the rock 
below the detector, and the contained flux when muons are created in the (ice) detector. We also 
calculate the event rates for showers from neutrino interactions in the detector and show that the 
signal dominates over the background for f 50GeV < m x < 1 TeV for = fOO GeV. 
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I. INTRODUCTION 

Dark matter's presence is inferred from gravitational 
effects on visible matter at astronomical scales. A wide 
range of observational data show that the dark matter is 
cold or warm (i.e. it became non-relativistic before or at 
the time of galaxy formation) and composes about 23% of 
the total density of the Universe [l[ . There are no viable 
candidates for dark matter within the standard model of 
elementary particles, but many in proposed extensions 
of the standard model. Among these, weakly interact- 
ing massive particles (WIMPs) of mass in the 100 GeV 
to several TeV range provide a natural explanation for 
the observed dark matter density @. We are going to 
concentrate on WIMPs in this paper. 

Although the detection of dark matter particles may 
be possible at the Large Hadron Collider (LHC), find- 
ing them in direct or indirect dark matter searches will 
be necessary to determine if they are indeed stable on 
cosmological timescales and how abundant they are at 
present [3] . Many direct or indirect dark matter searches 
are being carried on at present Q. Indirect dark mat- 
ter searches look for WIMP annihilation (or sometimes 
decay) products, either photons [5;-7] or anomalous cos- 
mic rays, such as positrons and antiprotons (8rll3|, or 
neutrinos 

EMI. 

For some years, observations of an 
excess in the positron fraction e + /(e + + e~) by HEAT 
(the High Energy Antimatter Telescope) |9|, a bright 511 
keV gamma-ray line from the Galactic Center by INTE- 
GRAL (the International Gamma Ray Astrophysics Lab- 
oratory) Q and a possible unaccounted-for component of 
the foreground of WMAP around the galactic center, the 
"WMAP Haze" [7] (among others) have been considered 
possible hints of WIMP dark matter annihilations. 

More recently, the PAMELA satellite (Payload of 
Antimatter Matter Exploration and Light-nuclei Astro- 
physics) reported an excess in the positron fraction in 



the energy range of 10-100 GeV with respect to what 
is expected from cosmic rays secondaries [T(| , which con- 
firmed the HEAT excess. Also ATIC (the Advanced Thin 
Ionization Calorimeter) and PPB-BETS (the Polar Pa- 
trol Balloon and Balloon borne Electron Telescope with 
Scintillating fibers) observed a bump in the e + +e~ flux 
from 200 to 800 GeV [HEl, but this was not confirmed 
by the air Cherenkov telescope HESS [l3j and by the 
Fermi Gamma Ray Telescope. Fermi found a slight ex- 
cess in the e++e~ flux between 200 GeV and 1 TeV 0. 

Indirect searches for dark matter annihilations via neu- 
trinos with experiments such as AMANDA (Antarctic 
Muon And Neutrino Detector Array) [l5[ and IceCube 
[l6| also constrain dark matter models. The cubic kilo- 
meter size neutrino telescope (KM3NeT), planned to be 
built at the bottom of the Mediterranean Sea pjj, will 
provide additional constraints, with its different view of 
the sky and in particular, the galactic center. Many the- 
oretical studies have concentrated on the indirect dark 
matter detection via neutrino signals 

The positron excess observed by PAMELA may be 
explained by the presence of particular astrophysical 
sources (g.S^ pulsars) [23| , or by the annihilation [24], [H| 
or decay [26| of dark matter particles. If the observed 
anomalies in the PAMELA and FERMI data are due 
to dark matter annihilation, a larger annihilation rate 
than expected for typical thermal relics must be assumed. 
This enhancement may happen due to either large inho- 
mogeneities in the dark matter distribution near Earth 
(subhaloes) and/or a larger annihilation cross section of 
the dark matter particles. This last possibility may hap- 
pen if the dark matter particles are not thermal relics 
P, HH, in which case they can have larger annihilation 
cross sections in the early Universe, or due to an en- 
hancement of the annihilation cross section only at very 
low velocities [27j . which would not affect their annihi- 
lation in the early Universe. Whatever its origin may 
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be, the needed enhancement is quantified by a "boost 
factor," B, ranging from 10 to 10 4 The typ- 

ical WIMP thermal relic annihilation cross section is 
(av) = 3 x 10~ 26 cmV 1 . 

WIMP models explaining the PAMELA positron ex- 
cess must be peculiar in other aspects as well. To avoid 
overproduction of antiprotons, the dark matter annihila- 
tion or decay must proceed dominantly to leptons. More- 
over, the absence of a sharp shoulder in the electron plus 
positron spectrum (that had been observed by ATIC) in 
the Fermi data corresponding to an energy close to the 
parent dark matter particle mass means that the direct 
production of electrons must be suppressed with respect 
to the production of electrons (and positrons) as secon- 
daries. Final states including r's or /i's of dark matter 
not lighter than 1 TeV fit the PAMELA, HESS and Fermi 
data best [28j . These leptophilic dark matter candidates 
f24| would copiously produce neutrinos [l9| whose fluxes 
are constrained by the observations of Super Kamiokande 
(SK) 2j| toward the direction of the Galactic Center. 
Neutrinos with energies of the order of the dark matter 
mass, E v < m x , would propagate without being deflected 
towards the Earth. However, during their travel, vacuum 
oscillation effects would mix the three flavors. Some frac- 
tion of the arriving muon neutrinos would be converted 
into muons via charged current interactions in the Earth 
which can be detected in Earth based neutrino telescopes. 

Neutrino signals in underground or underwater detec- 
tors of dark matter annihilation in the Galactic Center 
are the subject of this paper. We calculate the neutrino 
induced upward and contained muon flux, as well as the 
neutrino induced muon and shower event rates due to 
dark matter annihilation in the Galactic Center. We take 
into account the muon propagation in the Earth when 
evaluating the upward muon flux [30| and study the en- 
ergy range of muons for which upward muon events dom- 
inate over the contained ones. We show that the shape of 
upward muon fluxes differs significantly from the shape 
of the neutrino spectra at production, due to the smear- 
ing produced by neutrino interactions and muon propa- 
gation. The muon propagation shifts the flux to lower 
energies, while the contained muon flux increases with 
muon energy due to the linear energy dependence of the 
neutrino charged-current interaction. We consider differ- 
ent WIMPs annihilation channels that contribute to the 
neutrino signal, including direct annihilation to neutri- 
nos, to charged leptons and to quarks or gauge bosons. 
We evaluate rates of contained events and upward events, 
of relevance to IceCube and future neutrino detectors like 
KM3NeT. 

In the next section, we evaluate expressions for muon 
flux from the incident neutrino flux interacting with the 
medium. In Section III we present our results for muon 
flux and muon event rates from the annihilation of the 
dark matter in the Galactic Center compared with the 
atmospheric background and evaluate rates for hadronic 
and electromagnetic showers. Finally, in Section IV we 
summarize and discuss our results. 



II. MUON FLUX 

The neutrino flux at the Earth due to the annihilation 
of dark matter particles with mass m x in the Galactic 
Center is given by 

^=*xfeV^| (1) 
dE v \f£ F dE v J w 

where R is the annihilation rate given by: 

R = B^%fdnf dl(9)p 2 (l) 7 (2) 
8nm x J J l o s 

dN^ jdE v is the neutrino spectrum at the production for 
a given annihilation channel F with branching fraction 
Bp, B is the boost factor, p(l) is the dark matter den- 
sity, integral is over the line of sight (l.o.s) within a solid 
angle AO, centered in the Galactic Center. The neutrino 
energy distribution, dN v jdE v , depends on the particle 
produced. Some examples appear in Appendix A. For 
all of the evaluations below, we take the dark matter an- 
nihilation cross section to have the typical thermal relic 
value (av) = 3x 10~ 26 cm 3 s _1 . 

For practical reasons the dimensionless quantity (J2)si 
is defined in which the dark matter density profile p(l) is 
embedded [l9j . 
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where 1(9) is the distance from us in the direction of 
9 which is the cone half angle from the Galactic cen- 
ter, R is the distance of the solar system from the 
Galactic Center and p is the local density near the so- 
lar system, which are taken to be R = 8.5 kpc and 
p = 0.3 GeVcm -3 . As a practical matter, we consider 
two profiles, the Navarro-Frenk- White (NFW) 31] profile 
and a cored isothermal profile. Some typical values for 
(J2)nAO can be found in Ref. 132], where (J2)siAO= 
6.0(10.0) for 9 = 5° (10°) for the NFW profile, and 
(J 2 )nAQ= 1.3(4.3) with 9 = 5°(10°) for the isothermal 
profile. 

The high energy neutrinos coming from the Galactic 
Center then interact with the matter in the Earth and 
produce muons that traverse to the detector (upward 
events), or they interact in the detector producing muons 
or showers (contained events). Muon range or stopping 
distance, R^E^, E th ), is given by 

1 / a + BE], \ 

w e ->=s4^j <4 » 

where a corresponds to the ionization energy loss and 
B accounts for the bremsstrahlung, pair production and 
photonuclear interactions. For example, for a muon with 
initial energy E^ ~ 1 TeV, when E t h — 1 GeV the muon 
range is roughly 1 km whereas the decay length of a muon 
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with the same initial energy is much larger (~ a few thou- 
sand kilometers). For detectors with a characteristic size 
of 1 km 3 , contained events are most important for WIMP 
masses below about 1 TeV, while for smaller detectors 
like SuperK, upward events are relatively more impor- 
tant. 

Using Eq.(l) and following the theoretical framework 
presented in Ref. [30|, the upward muon flux at the de- 
tector is given by 
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Here P SU rv accounts for muon energy loss in transit from 
its production position to the muon's entry into the de- 
tector. For an energy independent energy loss parameter 
P, the survival probability is 



Psurvi-E^, E^) ~ I 



(6) 



where T = m^j '(cr^ap) in terms of the muon mass, muon 
lifetime and the density of the medium p in g/cm 3 . 

For production in the detector, the contained muon 
flux is 
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where D is the size of detector. The quantity dPcc is 
the probability for a neutrino with energy E v to convert 
into a muon within the energy interval of dE 1 ^ and over 



a distance dz: 
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where Na — 6.022 x 10 3 is Avogadro's number. The 
differential cross sections da^ n jdE 1 ^ are the weak scat- 
tering cross sections of (anti-)neutrinos on the nucleons, 
which can be approximated by [33| 



d(J l^ _ 2m p G 2 F 
dEi ~ 7T 
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(9) 



the parameters a and b for charged current scattering are 
shown in Table Q] 

Muon rates, A^(m x ), are obtained by integrating 
Eqs.(5) and (7) over the muon energies, i.e., 
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TABLE I: Parameters for the charged current neutrino- 
nucleon differential cross section, as noted in Ref. [331. 
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TABLE II: Parameters for the neutral current neutrino- 
nucleon differential cross section, as noted in Ref. 1331 - 



where E t h is the muon detector threshold. 

Another set of possible signals of dark matter are 
the showers produced in neutrino charged-current and 
neutral-current interaction in the detector. The con- 
tained shower flux in CC and NC interactions is given 

by |33: 
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where the shower energy is 
E s h ~ E v 
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The neutral current cross section can also be approxi- 
mated with Eq. © where the parameters a and b appear 
in Table HU 

In the limit of the survival probability P surv going to 
unity, the energy dependent flux can be calculated ana- 
lytically when Eq. (9) is used for the neutrino-nucleon 
cross section. The analytic results for a variety of decay 
channels are shown in Appendix B. 



III. RESULTS 

The direct production channel, \X ~^ v^v^, where \ is 
the WIMP, is the most promising channel for the detec- 
tion of dark matter annihilation, assuming an adequate 
annihilation cross section, because of the monoenergetic 
neutrinos. A typical example of a dark matter particle 
candidate which annihilates into a neutrino pair is the 
lightest Kaluza-Klein particle. However, some particle 
candidates, for example neutralinos and leptophilic dark 
matter, produce neutrinos only as secondary particles, 
via the decay of the particles into which the dark matter 
particles annihilate, such as p + \i~ , t + t~, 66, W + W~, 
etc. 

In the first two figures, we present our results for the 
differential upward muon flux due to the annihilation of 
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TABLE III: Parameters for the atmospheric and Vp flux, 
in units of GeV _1 km _2 yr _1 sr _1 . 



a dark matter particle via the direct production (xX ~~ ^ 
v iiPii) channel. To illustrate various contributions, we 
choose the dark matter particle mass m x =500 GeV, and 
for Fig. 1, the NFW dark matter density profile [3l| and 
the boost factor B =200 which is in the range of the 
boost factor values that explain the PAMELA data (2lj . 
For Fig. 2, the dark matter density profile is the cored 
isothermal profile and we use a boost factor B = 800 to 
match the normalization of the NFW density profile for 
the 5° cone half angle. 

We show our results for two different choices of the 
cone half angle (5° and 10°) and compare them with 
the angle-averaged background due to the atmospheric 
neutrinos (in units of GeV _1 km _2 yr _1 sr _1 ) 
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which was obtained using the angle- dependent atmo- 
spheric neutrino flux parametrization in Ref. [35] , 
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The values of the parameters No, 7, a, b, c and e, given 
in Table III, were determined by fitting angle-dependent 
atmospheric neutrino data from Ref. 3(i . The resulting 
final muon flux with this approximated neutrino back- 
ground is about 50% larger (smaller) than that from the 
vertical (horizontal) atmospheric neutrinos. 

For a 10° cone half angle, the signal dominates over 
the background in the range 180 GeV< E^ <420 GeV 
for the NFW profile. We note that the background signal 
is suppressed more than the dark matter signal with the 
decrease in the cone of half angle. As a comparison, for 
a 5° cone half angle the signal exceeds the background in 
a wider range of energies, 60 GeV< <480 GeV. 

From Fig. 2, we note that in case of the isothermal 
profile for the dark matter in which there is a relatively 
less dense core region in the isothermal profile, by in- 
creasing the cone half angle from 5° to 10°, there is an 
almost equal enhancement of the upward muon fluxes 
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FIG. 1: Upward muon flux obtained from dark matter an- 
nihilation into neutrinos in the Galactic Center, for a cone 
half angle (6) of 5° (dot-dashed) and 10° (dotted). The back- 
ground upward muon fluxes due to (angle-averaged) atmo- 
spheric neutrinos are shown with the solid (for 8 = 5°) and 
the dashed (for 9 = 10°) curves. The NFW dark matter pro- 



file is used, along with a boost factor B 
GeV. 



200 and m y 



500 



from the atmospheric neutrino background and from the 
dark matter annihilation in the center of the galaxy. For 
the set of the parameters that we choose here, the dark 
matter signal becomes larger than the background in 
the energy ranges of 100 GeV< E^ <470 GeV and 70 
GeV< Ef, <480 GeV for the cone half angles 10° and 5°, 
respectively. 

Fig. 3 shows the dependence of the differential muon 
fluxes from dark matter annihilation via the direct pro- 
duction channel for m x = 200, 500 and 800 GeV. We 
again consider the NFW profile, a fixed boost factor 
(B = 200) and a fixed cone half angle (6 = 5°). The 
figure shows the upward flux as well as the contained 
flux assuming a detector size D — 1 km in Eq. (7). 
We find that regardless of the mass dependence, the up- 
ward event spectrum is a decreasing function of the muon 
energy whereas the corresponding spectrum of the con- 
tained events increases with the muon energy up to the 
cut-off set by the initial neutrino energy. In our calcula- 
tions, we assume that the dark matter particles annihi- 
late at rest and thus the neutrino energy for this decay 
mode can be set to the rest mass of the dark matter 
particle, E^ = m x . 

The signal for the muon flux from the contained events 
has a stronger suppression with the increase in the dark 
matter mass than for the upward muon events. This is 
due to the m~ 2 dependence in Eq. (2). The mass de- 
pendence for upward events is more complex because of 
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FIG. 2: Same as Fig. 1 but for the cored isothermal dark 
matter density profile and a boost factor B = 800. 

the mass dependence in the upper limit of the z integra- 
tion in Eq. (5). A large mass m x (and therefore higher 
E v ) produces a higher energy muon which has a longer 
range in the rock below the detector. For example, for 
Eft > 380 GeV, the upward event signal from the anni- 
hilation of the dark matter particle with mass m x = 800 
GeV dominates over the one from that of the dark matter 
particle with mass m x = 500 GeV. 

For a wide range of muon energies, the dark matter 
signal is above the atmospheric background both for con- 
tained and upward events in the \X ~~ ^ v iiPii channel with 
the boost factor used here. We find that for a given dark 
matter mass the contained events exceed the upward ones 
in the range E^ > 0.6m x . 

In Fig. 4, we present our results for the differential 
muon flux due to XX — > t + t~ channel. This channel 
is characteristic of all three-body decays into neutrinos 
(secondary neutrinos) . Again shown are the upward and 
contained signals from m x = 200, 500 and 800 GeV with 
the NFW profile and B — 200. 

Note that in the case of secondary neutrinos, the sig- 
nal for both upward and contained events decrease as the 
muon energy increases, and for a fixed m x , the contained 
events, in general, dominate over the upward events for 
muon energies lOOGeV < E^ < m x . This is a conse- 
quence of considering a detector size of D = 1 km, a size 
larger than the range of a muon with an energy of less 
than 1 TeV. The figure shows that even for a half angle 
of 5°, in case of NFW profile one would need a boost 
factor on the order of about 2000 for the dark matter 
signals from the secondary neutrinos to be above the at- 
mospheric background. 

Measurement of the muon flux can also be used to 
distinguish different dark matter models, as seen in Fig. 



FIG. 3: Muon flux due to the dark matter annihilation into 
neutrinos in the Galactic Center for different dark matter 
masses, curves correspond to the dark matter masses of 200 
GeV, 500 GeV and 800 GeV, respectively. The corresponding 
backgrounds are also shown. All the solid lines correspond to 
the contained events with D = 1 km, whereas the dashed ones 
to upward events. 



5 where we compare signals from different annihilation 
channels: XX ~^ W + W~ , XX ~~ ^ t + t~ and XX bb for 
the NFW profile, with B = 200, the half angle equal 5° 
and m x = 500 GeV. The signals from the b-quark and 
the tau decay modes differ only by an overall factor which 
is close to the ratio of the decay branching fractions of the 
corresponding modes given in the Appendix I. However, 
for the W decay, being a 2-body decay, the shape of the 
differential muon spectrum is quite different than those of 
the b-quark and tau which are both 3-body decay modes. 
This indicates that muon flux from the secondary neu- 
trinos as a by-product of the dark matter annihilation 
can also be useful in discriminating different dark matter 
models. 

We now turn to the total rate of upward and contained 
muons produced by + from direct dark matter an- 
nihilation to neutrinos. Integrating the differential fluxes 
over the final muon energy, we obtain the muon rate 
from the annihilation of the dark matter as a function of 
the mass m x (Fig. 6) for the NFW profile with B = 200 
and 9 = 5°. Here, the threshold energy is taken to be 
Eth — 80 GeV. Due to the finite size of the detector 
[D = 1 km), and m~ 2 dependence of the annihilation 
rate, the signal for the contained events decreases with 
increasing the dark matter mass. On the other hand 
for upward events, heavier dark matter particles yield 
more energetic neutrinos which makes a larger portion of 
muons in the rock below the detector to contribute to the 
final muon flux. This effect combined with the energy de- 
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FIG. 4: Muon fluxes due to the secondary neutrinos produced 
through the dark matter annihilation into tau particles in 
the Galactic center for different dark matter masses; m x — 
200, 500 and 800 GeV. The solid (dashed) curves correspond 
to contained (upward) events. 



pendence of the neutrino charged-current cross section, 
results in increasing muon rate up to m x — 650 GeV, 
at which point the ni~ 2 dependence of the annihilation 
rate takes over resulting in slow decrease of the muon 
rate. Comparison of contained and upward muon rates 
presented in Fig. 6 indicates that for m x < 500 GeV 
the signal from the contained events still dominates over 
the signal from the upward events. Even though the sig- 
nal depends weakly on the value of the threshold energy, 
the background is very sensitive to it due to the large 
contribution from the low energy atmospheric neutrinos. 
The signal to background ratio increases with increasing 
the muon energy threshold. We obtain the same results 
for the isothermal dark matter density halo profile if the 
boost factor is taken to be 800 for the same cone half 
angle of 5°. 

In Fig. 7 we show our results for the 10° cone half 
angle. We note that in case of contained events the signal 
dominates over the background for 100 GeV < m x < 200 
GeV, when the threshold energy is 80 GeV. For upward 
events, signal is below the background for all m x . The 
isothermal dark matter density halo profile gives larger 
signal than obtained with the NFW profile by about a 
factor of 2, due to its larger increase of (J2)n for 10° 
relative to 5°. 

In Fig. 8 we show contour plots for upward muon 
events, = (0.5, 5, 50, 500, 850)km^V _1 - The solid 
(dashed) lines correspond to the muon energy thresh- 
old of 50 (80) GeV. We also calculate that iV M = 
714(516)km _2 yr _1 for the upward muon events due to 
the atmospheric muon neutrinos for the muon energy 



FIG. 5: Muon fluxes due to the secondary neutrinos produced 
through the dark matter annihilation into W bosons, tau par- 
ticles, and bottom quarks in the Galactic Center. The solid 
(dashed) lines for each channel correspond to contained (up- 
ward) events. The detector size is taken to be D — 1 km, 
and the cone half angle is 6 — 5° for the NFW profile with 
B = 200. 



threshold of 50 (80) GeV. We find that for a fixed cone 
half angle the annihilation cross section does not depend 
on m x for m x > 200 GeV to produce a given total muon 
flux since the decrease in the annihilation rate with m x 
is compensated with the increase in the muon range and 
neutrino cross section with m x . The dependence on the 
choice of the threshold is also negligible. However, for 
low mass dark matter particles, higher values of the an- 
nihilation cross sections are required in order to have the 
same total muon flux. This is due to the fact that the 
neutrinos originated from this low mass dark matter an- 
nihilation mostly contribute to the muon flux at ener- 
gies less than the thresholds we choose. The parameter 
space above the dotted line is excluded at 90% C.L. by 
Super-Kamiokande observations toward the direction of 
the Galactic Center with a cone half angle of 5° [29]. 

The dominant atmospheric neutrino flavor at neutrino 
energies above 40 GeV is Vy, which produces track-like 
events through charged current interactions in the neu- 
trino telescopes. Identifying track- like events could re- 
duce the background substantially. Recently it has been 
argued that IceCube+DeepCore will be able to put con- 
straints on dark matter properties in a more efficient way 
by just analyzing the cascade (i.e shower) events which 
are due to charged current interactions of v e%r and the 
neutral current interactions of the all neutrino flavors 
[22| . Since the weak scattering cross sections are inde- 
pendent of the flavors, the signal-to-background ratio is 
enhanced in shower events since can only contribute 
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FIG. 6: Total muon fluxes due to the dark matter annihila- 
tion into neutrinos in the Galactic Center. The solid (dashed) 
lines correspond to contained (upward) events. 



FIG. 8: Upward muon events curves, = 
(0.5, 5,50,500, 850)km" 2 yr"\ for the energy threshold 
of 50 GeV and 80 GeV are shown by the solid and the dashed 
lines, respectively. The boost factor is set to be unity and 
the cone half angle is chosen to be 5°. 
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FIG. 7: Same as Fig. (6) but for 10° cone half angle. 



to the shower events through neutral current interactions 
where the cross section about 1 /3 of the charged current 
cross section. 

In Fig. 9, we show hadronic shower rates as a function 
of m x from neutral-current and charged-current interac- 
tions of muon neutrinos and antineutrinos. These rates 
are the same for any other neutrino flavor with a demo- 
cratic XX ~~ ^ y v annihilation rate. Also shown is the 
hadronic shower rate due to the atmospheric muon neu- 
trinos; JV s ^ m = 524(168)km~ 2 yr~ 1 for the charged cur- 



rent (neutral current) interactions. The shower threshold 
is taken to be 100 GeV. We note that the background due 
to the atmospheric electron and tau neutrinos is much 
smaller than for the muon neutrinos, so the signal to 
background would not change much here when all the 
neutrino flavors were included. 

We also evaluate the electromagnetic shower rate as 
a function of m x due to electrons produced by the 
charged-current interactions of v e , with an electromag- 
netic shower threshold set at 100 GeV. The atmospheric 
shower rate is evaluated using the atmospheric v e and I7 e 
flux for an effective zenith angle 0.4 < cos 9 Z < 0.5, which 
roughly corresponds to the angle describing the position 
of the Galactic Center relative to the IceCube, 



d<f> \ _ 500.0 / E 

dEdil J Ve ~ (GeVm 2 ssr) \ GeV / 

d<p \ _ 382.6 / E 

dEdQ )- e ~ (GeVm 2 ssr) I GeV 



-3.57 



-3.57 



(15) 



From Fig. 10 we see that the signal-to-background ra- 
tio is increased for the electromagnetic showers relative 
to hadronic showers (see Fig. 9) mainly due to a very 
small atmospheric electron neutrino flux which is about 
34km~ 2 yr~ 1 . For secondary electron neutrinos from the 
decay of taus which are produced via XX ~^ t + t~, the 
signal becomes comparable to the background. 

For the future neutrino detector which is positioned 
in the northern hemisphere, such as KM3Net, the rele- 
vant background would be coming from almost horizon- 
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tal showers, which is about a factor of three to four times 
larger than the flux given by Eq. (14), giving approxi- 
mately electromagnetic shower flux of 100 km - yr _1 . 
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FIG. 9: Hadronic shower rates for charged-current (dashed) 
and neutral current (dot-dashed) interactions of u^+V^ when 
muon neutrinos are produced directly from the dark matter 
annihilation in the Galactic Center, compared with the atmo- 
spheric background. The NFW profile, with B = 200, 6 = 5° 
and D — 1 km are used. 

In Fig. 11 and 12, we present the contour plots for 
contained showers with the energy threshold of 100 GeV. 
The main difference between the showers and the upward 
muons appears for m x > 200 GeV where for a given total 
number of shower events higher annihilation cross sec- 
tions is required with the increase in m x . This is due to 
the contained event nature of the shower events which are 
all produced inside the detector with finite size. Thus, in 
contrast to the case for the upward muon events that we 
discussed earlier, the strong suppression of the annihila- 
tion rate with m x can not be compensated because of the 
finite size of the detector. The charged current showers 
actually require a smaller annihilation cross sections in 
order to produce the same number of total shower events 
that neutral current showers produce for a fixed m x due 
to the larger weak scattering cross sections. 

The signal detection significance can be evaluated us- 
ing 



S = 



V s 



y/(N s + N b ) ' 



(16) 



where N s corresponds to the number of events for the 
signal, while Nb is the background. We obtain the time 
it would take to observe a 5a effect using our results for 
the contained muon events (Fig. 6), hadronic showers 
(Fig. 9) and electromagnetic showers (Fig 10), 
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FIG. 10: Electromagnetic shower rates as a function of m x for 
v e + v e charged-current interactions when electron neutrinos 
are produced directly in the annihilation of dark matter in the 
Galactic Center, compared with the atmospheric background 
for shower energies above 100 GeV. The NFW profile, with 
B = 200, 6 = 5° and D — I km are used. 



t 



25(V S + N b ) 
N?V 



(17) 



where V = 0.04(0.02)km 3 is the effective volume of Ice- 
Cube+DeepCore for the track-like (shower) events. In 
Fig. 13, we show the observation time (t) required for Ice- 
Cube+DeepCore detector to detect or exclude the dark 
matter signal via the direct production channel at a 5a 
level. Here, we again use fixed boost factor (B = 200) 
and cone half angle (6 = 5°). Our results, when we take 
Bf = 1 for the direct production channel, suggest that 
in less than two years of observation IceCube+DeepCore 
will be able to reach a 5a detection for the contained 
muon and electromagnetic shower events for a wide range 
of m x . Decreasing the branching fraction by an order 
of magnitude increases the observation time significantly 
in order to reach the same significance. For instance, 
t ~ 10 - 50 years, for 150GeV > m x < 500GeV in the 
case of contained muon events, and somewhat shorter for 
the electromagnetic showers. 

In the case of secondary neutrino production, when 
neutrinos are produced from tau decays, and taus are 
products of dark matter annihilation, these neutrinos can 
interact inside the detector producing hadronic and elec- 
tromagnetic showers, in addition to muon neutrinos pro- 
ducing muons via charge-current interactions. In Fig. 
14 we show that IceCube+DeepCore detector could po- 



tentially detect a 2a effect in 5 (8) years for m x 



300 



GeV (ITeV), in case of excluding muon-like events. To 
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FIG. 11: Hadronic shower events curves, N^ h — 
(0.5, 5, 50, 500)km _2 yr _1 , for the charged current (solid) and 
neutral current (dashed) processes, for a NFW dark matter 
density profile, a 5° cone half angle, the boost factor set to 
be unity and D — 1 km. 

reach a 2a detection for the electromagnetic showers due 
to the secondary electron neutrinos IceCube+DeepCore 
will need about 10 — 20 years of observation for 250 GeV 
< m x < 1 TeV. When muon-like events are included, 
the observation times for the hadronic showers become 
similar to those for the electromagnetic showers. The 
time needed for a 5a effect for hadronic (electromagnetic) 
showers is almost an order of magnitude longer than for 
a 2cr effect. 

Comparing the secondary and direct production (Fig. 
13) one sees that it takes longer (by about one order of 
magnitude) to detect showers from secondary neutrinos 
that to detect showers from primary neutrinos. This is 
because of the different shape of the shower energy dis- 
tributions: for direct neutrinos it increases with energy 
and for secondary neutrinos it decreases with energy. 

Since the angular resolution for showers is expected to 
be much worse than for muons, for the angular resolu- 
tion of 30°, the number of signal events will be larger 
by a factor of 6, while the background will increase by a 
factor of 35, which results in reducing the time it would 
take IceCube+DeepCore to see a 2a effect to 3 years for 
hadronic showers without track-like events. This is in 
qualitative agreement with the results presented in Ref. 

For dark matter models in which neutrinos are decay 
products of taus produced in the dark matter annihi- 
lation, looking for contained hadronic showers in Ice- 
Cube+DeepCore seems promising to detect a signal at 
the 2 sigma level, assuming the NFW dark matter halo 
profile and a boost factor B = 200. 



FIG. 12: Electromagnetic shower events curves, N*™ = 
(0.5, 5, 50, 500)km" 2 yr _1 for a NFW dark matter density pro- 
file, a 5° cone half angle, the boost factor set to be unity and 
D — 1 km. 

In Table (|IV[) we give a summary of our results for 
the event rates for various dark matter masses. We 
consider the direct production of neutrinos (xx ~* v v) 
and the neutrinos from the tau decay (xx ~* t + t~ 
l + l~v T v T vii>i). We classify the event rates as contained 
(ct) and upward (up) for the track-like muon (/i) events, 
and depending on the type of the interaction involved 
charged current (CC), neutral current (NC) and electro- 
magnetic (em) for the shower events. Two different cone 
half angles are chosen, 9 — 5° and 8 = 10°, and the 
threshold energy for the track-like muon (shower) events 
are set to be 80 (100) GeV. We also show the atmospheric 
neutrino background for the track-like muon and for the 
shower events. 



IV. CONCLUSION 

We have studied neutrino signals from dark matter an- 
nihilation in the Galactic Center. We have calculated 
contained and upward muon fluxes from neutrino inter- 
actions, when neutrinos are produced in annihilation of 
dark matter either directly or via the decay of taus, W- 
bosons or b-quarks. We have shown that in the case of 
direct neutrino production, the signal is above the at- 
mospheric background for both contained and upward 
events, assuming that the annihilation rate is enhanced 
by boost factor of 200 (when the NFW dark matter halo 
profile is used) and that the branching ratio of dark mater 
annihilation into neutrinos is one. In general, the boost 
factor values that are required to explain the data ob- 
tained by the indirect detection experiments vary de- 
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FIG. 13: Time as a function of dark matter mass, m x , for 
the direct neutrino production channel (xx to reach 

a 5ct detection level for IceCube+DeepCore detector. The 
curves correspond to hadronic showers (solid for neutral cur- 
rent, dashed for charged current interactions), electromag- 
netic showers (dotted) and the contained muon events (dot- 
dashed). Bf = 1(0.1) for the lower (upper) curves, the boost 
factor is taken to be 200 and the cone half angle is 5° for all 
curves. 



pending on the dark matter model and the dark matter 
mass. For the specific dark matter model our results can 
be rescaled by the corresponding product of the boost 
factor B and the branching ratio Bf- 

We have found that the contained muon flux domi- 
nates over the upward muon flux for all energies when 
m x = 200 GeV. However, as we increase the mass m x 
of the dark matter particle, for example when m x = 500 
GeV, the upward muon flux dominates up to = 300 
GeV, and for m x = 800 GeV, up to = 500 GeV. This 
is due to the increasing muon range as the muon ini- 
tial energy increases, which becomes possible when m x 
is larger thus producing higher energy neutrinos in the 
annihilation. In the case of secondary neutrino produc- 
tion, the signal becomes comparable to the background 
if the boost factor is an order of magnitude larger than 
the value we considered. We have shown that the shape 
of the muon flux depends on the specific decay mode, 
and that the dominant flux comes from tau decay at 
low muon energies, and from W-decay for muon energies 
above 200 GeV. The total upward muon rates have a 
weak dependence on m x and on the muon energy thresh- 
old for m x > 400 GeV, due to the balance of the energy 
dependence of the muon range, the upper limit of the 
muon energy (given by m x ) and the explicit dependence 
on m x (~ m~ 2 ) of the muon flux. However, the total 
contained muon rates show a sharp decrease with m x for 



FIG. 14: Time as a function of dark matter mass, m x , for 
the secondary neutrino production channel XX ~^ t + t~ — > 
l + l~v T P T i>ii>i to reach a 2a (solid curves) or a 5a (dashed 
curves) detection level when measuring electromagnetic show- 
ers (top curves) and hadronic showers without charged track- 
like events (lower curves). 



m x > 150 GeV due to the finite size of the detector. Up- 
ward muon events dominate over contained muon events 
for m x > 550 GeV. 

We have also shown that showers produced by neu- 
trino interactions, when neutrinos are produced directly 
in dark matter annihilation, could also be used to de- 
tect a dark matter signal from the Galactic Center. In 
particular, electromagnetic showers have much smaller 
background, from atmospheric electron neutrinos, than 
the hadronic showers. In addition, we have studied the 
contour plots of both the upward muon events and the 
showers and we have shown the required dependence of 
the annihilation cross section on the dark matter mass in 
order to observe a fixed number of event rates. We have 
discussed the origin of different shapes for the contour- 
curves in each case and pointed out the contained event 
nature of the shower events. We have shown that after 
one year IceCube+DeepCore detector could potentially 
observe a 5cr signal effect by measuring contained muons 
(for direct neutrino production), or in 5 to 8 years a 2a 
effect with hadronic showers even in the case when they 
are due to secondary neutrinos. 

IceCube+DeepCore will be able to identify track-like 
events due to the charged current interactions of muon 
neutrinos, the showers due to neutral current interac- 
tions of all the neutrino flavors and the charged current 
interactions of electron and tau neutrinos. In particular, 
above the neutrino energy of 40 GeV the signal to back- 
ground ratio for showers is further enhanced since the 
atmospheric tau and electron neutrino fluxes are sup- 
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TABLE IV: Event rates per km 2 per yr for the contained 
(ct), upward (u) muons (/x) and for the showers (sh) pro- 
duced via charged current (CC), neutral current (NC) and 
electromagnetic (em) interactions. Neutrinos from direct pro- 
duction (xx ~* v ~v) channel and secondary neutrinos from 
XX t + t~ channel are considered. We have set B-Bf = 200 
for each channel. The cone half angle is chosen to be 5° and 
10°. The threshold energy for the muon (shower) events is 
set to be 80 (100) GeV. The backgrounds due to atmospheric 
neutrinos are also presented. 



good angular resolution [2fj| . 

Due to its location in the northern hemisphere, the fu- 
ture KM3NeT experiment will be complementary to Ice- 
Cube+DeepCore in searching for neutrino signals from 
dark matter annihilation in the Galactic Center through 
the observation of upward muon events. The atmospheric 
muon background at the KM3NeT will be suppressed 
significantly since the Earth will act as a shield to those 
muons. Independent searches of the upward muon events 
by KM3NeT and the contained muon and shower events 
by IceCube+DeepCore look promising for the discovery 
of the mysterious dark matter particle or for setting strin- 
gent constraints on its properties. 
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Appendix A: Neutrino Energy Distributions 

1. Neutrino energy distribution from direct 
production 

The neutrino energy distribution when neutrinos are 
produced directly from dark matter annihilation is given 
by a delta function, 



pressed relative to the atmospheric muon neutrino flux. 
Thus, the main background is the neutral current inter- 
action whose cross section is about a factor of three less 
than the charged current cross section of the atmospheric 
muon neutrinos. The measurement of the ratio of track- 
like muon and shower events eliminates the dependence 
on some parameters of the theory (e.g., boost factor, the 
dark matter density profile, etc) which only determine 
the overall normalization for the energy dependent dif- 
ferential muon fluxes, so the physical properties of the 
dark matter particle can better be determined. 

In addition to the boost factor due to Sommerfeld en- 
hancement that we have considered, there is potential 
enhancement of the dark matter signal due to the exis- 
tence of small substructures in the Milky Way Halo [38[ . 
Possible observation of this additional boost may be dif- 
ficult to observe because of the small population of these 
substructures unless the neutrino detectors have a very 



= 6(E„ - m x ) (Al) 

where the assumption is that the dark matter particles 
are essentially at rest when they annihilate. 



2. Neutrino energy distribution from r + r and bb 
decay modes 

In these decay modes, we use the unpolarized decay 
distributions, so the v and v distributions are assumed 
to be the same. The decay branching fraction is de- 
noted by Bf for a given decay mode /, / = r, b. The b 
quarks hadronize before they decay into neutrinos. The 
hadronization effect is taken into account by scaling the 
initial quark energy, Ei n = m Y , in the form Ef = Zftn x , 
where zt = 0.73 for b quarks [39j|. 
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The neutrino energy distribution from the decay of / = 
t + , t~ , b or b from \X ff i g approximately 

dN,, IB f . o o. , E,, 

^-(1 -3x 2 + 2x 3 ), where x = < 1 , 



dE„ E f 



Er 



(A2) 

where for each neutrino or antincutrino flavor (u e , v ei v^, 

f (m x , 0.18) r decay, 

J J | (0.73m x , 0.103) b decay . 

The energy distribution of the tau neutrinos from the 
decay of / = b or b is given by (1A2[) and the distribution 
from the decay of t + or t~ is given by [34| , 



dN„ 
dE v 



3Ef 



(1 — x 3 ), where x 



E± 
E f 



< 1 



(A4) 



3. W + W decay mode 

In the W + W~ mode, when the dark matter particle 
is at rest when it annihilates, Ew = m x /2 and ftw — 

1 — m 2 v /m x . The decay distribution, for each W, is 
dh v m x pw 2 2 

(A5) 

Here, -B = 0.105 for each neutrino flavor. 



Appendix B: Muon energy distribution 



The differential muon flux for the XX 
can be given as 



dE„ 



a/p) [ 



a(m x - E^) 



3m x 



K - E l) 



where 



B 



RoplBp (av) F (J 2 )nAnm p G 2 F N A 

4tt 2 p 



(Bl) 



(B2) 



There is a separate distribution for neutrino and antineu- 
trinos, since the parameters a and b depend on the inci- 
dent particle and the target. Here, for isoscalar nucleon 
targets, a = a v>v = 0.20,0.05 and b = b u<v = 0.05,0.20. 
Also appearing are the Fermi constant Gf — 1.17 x 10 -5 
GeV~ 2 and Avogadro's number Na — 6 x 10 23 . For stan- 
dard rock, a ~ 2 x 10 -3 GeVcm 2 /g accounts for the ion- 
ization energy loss and /3 ~ 3.0 x 10~ 6 cm 2 /g accounts 
for the bremsstrahlung, pair production and photonu- 
clear interactions and we take p = 2.6 g/cm 3 . 



For the contained events, a similar expression can be 
derived as 



dE^ m x 



--^(« + ^J©(^-^) (B3) 



where Q(x) = 1 if x > and Q(x) = otherwise, and 

c , = DB RqpIBf (av) F (J2)nAnm p G F N A p 
4tt 2 

where D is the size of the detector. 
We note that 



oc p° for the upward events 



dE, 



dE. 



- 1 for the contained events, 



oc p 



so, the muon flux doesn't depend on the rock density for 
the upward events except through a and /?, whereas for 
the contained events, the muon flux is directly propor- 
tional to the density of the medium. 

All the expressions for the muon flux derived below 
contain a <d(m x — E^) function. For secondary neutrinos 
which possess an energy spectrum in the form 



dN 
~dE 



A 



E„ 



(B5) 



where A is an overall factor, the differential upward muon 
flux can be calculated by using 



cA 



dE u 



(n+2) 



(E, + f ) 



-[P(m x ,£^,n) + K(m x ,E^,n) + 



+ L(m x ,E t _ l ,n) + M(m x ,E f _ l ,n)} 



(B6) 



vv channel where 



P(m x ,E fi ,n) 
K(m x ,E^n) 
L(m x ,£^,n) 



> +1) (m*-^) 



a(m x n ^> - Eir z >) 



Mim^E^ 



n = -- 



(n+1) 

(n+2) p (n+2) 

(n + l)(n + 2) 
bm x n ~ 1 \m x — E 3 ) 
3(n- 1) 
b(m ( x n+2) - 4" +2) ) 



(n- l)(n + 2) 



(B7) 



for n ^ 1 and when n = 1, 



d£7„ 3m 3 (£ M + f) 

„ / &\ 3aE„ml 
x \m x [a+ ] ' x 



+ E'l [ bin 



2 

a 6 

2 ~ 3 



(B8) 
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For the contained events and when ti^l, 



c'A r a 



- C ' A r fl ( m ("+l) _£(»+!))+ ^ ^ 2 * " 



dE» TO i"+ 2 ) L (n + l) v x 



bE 2 



+ -r-^irnp-V-EjT-V)] (B9) 

when 



(BIO) 



which reduces to 
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